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We consider the bremsstrahlung and model-dependent contributions to the radiative decay r~ — + 
, tt~tt v t j in the context of a meson dominance model. We focus on several observables related to 

f*"^ ■ this decay, including the branching ratio and the photon and di-pion spectra. Particular attention 

is paid to the sensitivity of different observables upon the effects of model-dependent contributions 
| and of the magnetic dipole moment of the p~(770) vector meson. Important numerical differences 

. i are found with respect to results obtained in the framework of chiral perturbation theory. 
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IT") , 1. INTRODUCTION 

Semileptonic tau lepton decays are a rich source of information about the properties of hadronic resonances below 
the tau lepton mass scale. They provide a clean environment to study the properties of charged /d(770) and ai(1260) 
i resonances which otherwise would be produced only through purely hadronic processes. The interplay of strong, weak 
and electromagnetic interactions in such processes offers an interesting place to test models for these interactions at 
low energies and to extract information about fundamental parameters of the standard model Q]. 

In this paper we are interested in the study of the radiative r~ — * Tr^7r°v T r y decay, a process that involves simul- 
taneously the three fundamental interactions at the lowest order. This decay channel has been studied previously in 
| references @>Q within different models and with different purposes. As is well known, the corresponding non-radiative 
42 ' T ~ ~~ * 7T _ 7r°zy T decay is dominated by the production of the p~(770) vector meson; thus, the emission of a single photon 
from this process is expected to carry information about the p--meson magnetic dipole moment 2). A meaningful 
v~( ■ extraction of this property from data is possible only with a full account of the model-dependent contributions to the 
radiative decay, which was not included in ref. 0- In this paper we pursue this study and consider the complete 
calculation of the radiative amplitude using a phenomenological model that includes all possible intermediate hadronic 
states. 

A different approach is followed in ref. |3( , where the radiative amplitudes were calculated in the framework of chiral 
perturbation theory and including resonances in the relevant kinematical regions. The interest of Ref. Q was focused 
on the relationship between the di-pion tau decay data and the its leading hadronic contribution to the anomalous 
magnetic moment of the muon a M jj| . As is known, present experimental information onr-> -k-kv decays are photon 
inclusive measurements Q. Thus, removing radiative effects from the measured di-pion mass distribution in such 
decays is important to predict the leading order hadronic vacuum polarization contribution to a M . A comparison of 
the two-pion mode in tau decays and e + e~ annihilations provides a sensitive test of the CVC hypothesis. At present, 
the prediction of a^ ad based on r — > nnv data seems to exceed by more than two standard deviations the corresponding 
prediction based on e + e~ data Q, even after the known sources of isospin breaking corrections are removed 
Since the production of high energy photons in t~ — ► n^ifiv^i decays is driven by the model-dependent contributions, 
a good account of the model-dependent effects is again mandatory. 



* Electronic address: 



afflores 


Eifis. cinvestav. mx 


glopezC 


8fis.cinvestav.mx 
fisica. unam.mx 


tolcdod 



2 



This paper is organized as follows: in section 2 we describe the necessary one-loop modifications of the propagator 
and electromagnetic vertex of the unstable p~ vector meson to achieve a gauge-invariant amplitude for the model- 
independent contributions; in section 3 we describe the form of the amplitude for the non-radiative r lepton decay and 
fix the parameters involved in our approximation; in section 4 we focus on the different contributions to the radiative 
decay amplitude and check their gauge invariance requirements; in section 5 we fix the coupling constants involved 
in the model-independent contributions and compute the different observables associated to the radiative two-pion 
t lepton decays; our conclusions are summarized in section 6 and an Appendix is devoted to discuss the kinematics 
associated to the this four-body decay. 



2. GAUGE INVARIANCE AND UNSTABLE PARTICLES IN RADIATIVE PROCESSES 

The physical amplitudes of radiative processes (Ai = e^Ai^, e being the photon polarization four- vector) have 
to satisfy the electromagnetic gauge invariance condition k^AA p — 0, where k is the photon four-momentum. As it 
has been discussed elsewhere jjj, when charged unstable particles are produced as intermediate states of a physical 
process some care must be taken to make compatible the unstable character of the resonance with the gauge invariance 
condition. One of the proposals to deal with this problem is the so-called fermion loop-scheme (fls) for gauge bosons 
||. According to the fls, only the fermion contributions in loop corrections to the propagator and electromagnetic 
vertex of gauge bosons have to be included to render gauge- invariant the resonant amplitudes || . 

In the case of hadronic resonances such as the p~ meson, it has been suggested that an analogous boson loop-scheme 
(bis) can be used to avoid such potential gauge pathologies. It has been shown 0,0 that when particles in loop 
corrections are massless, the corresponding dressed Green functions obtained in the fls and bis are the same as the 
ones obtained using the complex-mass prescription Mg — » M 2 — iMT (M and V are the mass and decay width of the 
resonance) in the bare Green functions. This prescription have been successfully used to describe experimental 
data of the elastic and radiative n + p scattering to extract the mass, width and magnetic moment of the A ++ baryon 
resonance. 

According to the boson loop-scheme, one has to include the absorptive parts of the one-loop corrections to the 
electromagnetic vertex and the propagator of the resonance in order to satisfy electromagnetic gauge invariance 0. 
In the case of a p~ vector-meson of mass m p and four-momentum q, the one-loop absorptive corrections arising from 
7r~7T° meson loops 1 gives the resonant propagator Q: 



D^_{q) = -i V VQ ' , (1) 

q 2 - jn 2 + i^/q 2 T p (q 2 ) 

where we have defined the energy-dependent width (in the limit of isospin symmetry m n - — m,o = m^) as follows: 

rp(9 2 ) = ^(? 2 -4m^( S -4m^), 

with gp-mr the pirn coupling constant (its value is discussed below). 

The one-loop absorptive corrections to the electromagnetic vertex (using the convention p~ a (p) — > p~^ (p')"f S (k) for 
Lorentz indices and four-momenta) gives the following result 0]: 

ieT al35 = ie{Tf s + Tf s ), (2) 

where 

rf 5 = {p + pT(f S + {k 9 aS - k 5 g afj )(3(0) ~ P P 9 a& ~ P' S 9 aP , (3) 

is the electromagnetic vertex at the tree-level, and /3(0) the value of the magnetic dipole moment of the p~ (770) 
meson in units of e/2m p (/3(0) = 2 corresponds to the normal or canonical value of the magnetic dipole moment; 
typical values of /3(0) computed in quark models lies in the interval 1.8 < /3(0) < 3.0 ^3)- The absorptive part of the 



1 The contribution of loops with K K° mesons can be included in a similar way, but we neglect its small contribution in this paper. 
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7r 7T° one-loop correction to the electromagnetic vertex of the p has been computed in ref Q using cutting rules. 
Its explicit form in the limit of isospin symmetry is given by: 



16ir(p 2 — p 12 ) 



A(p 2 )p a T 0S (p) - A(p' 2 )p' a T^(p') + B(p 2 )F af} (p)k s + B{p ,2 )F aS (p')k 



A{p 2 ) + B{p 2 )] \F af} (p)F^ {p) + F aS {p)F^{p) 



Pn 



where : 



A{p' 2 ) + B{p 12 ) F Q/ V)i^V) + F aS {p')F^(p') 

q 2 + (q 4 - Amlq 2 ) 1 / 2 



(4) 



B(q) = 2m^ln 

A(q 2 ) 
F^{q) = g^ 



Q 2 ~ (<7 4 — 4m 2 q 2 ) 1 / 2 
2(q 4 - 4 m \q 2 fl 2 



~{q 4 -±mlq 2 fl 2 



3q 4 
q^k v 



T^{q) = g^ 



qVq v 



(5) 



q-k q 

Since the above Green functions satisfy the Ward identity k a T a ps = [iDpsip)}^ 1 — [iDp${p')]~ l , the radiative 
amplitudes involving such vertices and propagators is automatically gauge-invariant. We will use this prescription in 
computing the radiative amplitude of t~ — > 7r~7r c V r 7 decay; as it will be discussed below, the model-independent 
contribution to this process involves the production and decay of an intermediate p _ (770) vector meson. 



3. NON-RADIATIVE TWO-PION DECAY 

In this section we focus on the meson dominance model for the non-radiative t~(P) — ► ir~ (Q)vr (Q')i/ T (P') decay, 
where the particles four-momenta are indicated within parenthesis. Our phcnomcnological model is based on the 
quantum-mechanical requirement of unitarity, according to which all possible intermediate states that are allowed to 
contribute given their quantum numbers have to be included (see Figure 1). In practice, only a few low-lying meson 
states are sufficient to describe experimental data. As it can be verified below, this model reproduces the Kiihn and 
Santamaria [llj parametrization of the vector form factor which contains the sum of the p~ (770) and of its higher 
excitations. 

In the limit of the isospin symmetry, the amplitude for this decay can be written in terms of a single vector form 
factor: 

M = ^^nQ-Q'),U(t) , (6) 

where Gp is the Fermi coupling constant, ^ M = u(P')^^(l — j^)u(P) denotes the leptonic current, t — (Q + Q') 2 is 
the square of the momentum transfer and V u <i is the Cabibbo-Kobayashi-Maskawa mixing matrix element. 

For the purposes of illustrating how the model works, we will assume that the amplitude is dominated by the 
exchange of two intermediate resonances: the p~ (770) and the p' (1450) vector mesons as shown in Figure 1. Applying 
the Feynman rules to the diagram of Figure 1 and using the vector-meson propagator given in eq. (1) we can obtain 
the following expression for the form factor: 



/+(*) 



m 2 -t~i VtTp (t) rn 2 , - t - my I> 



l + a [ m 2 _ l _ iVtT p (t) m 2 pl - i- im p >T p . 



(7) 



where g p (<y ) denotes the weak coupling of the p(770) (p'(1450)) vector meson (we neglect the energy dependence of 
the decay width of the p' meson) . 

The expression for the form factor in the second line of eq. (7), which coincides with the model of ref. [llj . 
follows from imposing the normalization condition f+(t = 0) = v2 and from the definition of the parameter a = 
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{m 2 p g p' g p' kit) / { m2 p i 9 p9 ptttt) ■ Using the experimental data on the p and p' decays 0] ( we take T(p' — > e + e ) = 1.48 
keV, T(p' — > 7r + 7r _ ) = 26.9 MeV, which we have estimated from relevant inputs in ref. ^3); we can obtain the 
following estimate for the relative strengths of these vector meson contributions in eq. (7): 



m p iT(p' — > e+e )T(p' — + im)(m 2 p — 4to^) 3 / 2 
y m p T{p — ► e+e~)r(;0 — * irTr)(m 2 , — Am 2 ) 3 / 2 
ps 0.102 . (8) 

This estimate is very close to the experimental value (|cr e:rp | = 0.120 ± 0.008) reported in ref. ^j|(see al so our n ^ 
discussed after eq. 10). This agreement renders confidence on the meson dominance model for the radiative decays 
to be discussed in this paper, and in particular about the values of the coupling constants extracted from other 
independent measured processes (see section 5). 

In order to provide a comparison with the results obtained for radiative r lepton decays in chiral perturbation 
theory hereafter we will restrict to the model with a single resonance, namely the p(770). In order to fix the 
parameters of this model, we have fitted the data of ref. [lj| for the pion form factor below Vt ps 1.1 GeV and have 
found that the approximation of using only one resonance gives a good description of data with the following central 
values for the resonance parameters: 

ra p = 776.66 MeV, g p ^ = 5.488 . (9) 
Using these values for the resonance parameters, we obtain the following result for the non-radiative branching fraction: 



B{t~ -> tTttV) = 20.75% (10) 

Clearly, our simple model underestimates the experimental value whose present world average is B exp (r~ — » ir~ir°is) — 
(25.47 ±0.13)% 0. This discrepancy can be attributed mainly to the fact that we have neglected the contribution 
of the p(1450) resonance which affects the higher energy tail of the hadronic spectrum (indeed, if we repeat the fit to 
data of reference using the two vector resonance model in eq. 7, fixing the mass and width of the p' to their PDG 
values 01 an( i assuming that a is real, we get m p = 775.80 MeV, g p7r7r = 5.867 and a = —0.12; this in turn leads to 
B(t~ — > ■n^-n a v T ) — 23.27%, which is closer to the experimental value). 

It is worth to mention that in ref. [3j, the authors have used the following form factor: 

f? P ® - ml _^ pT a P{i) *M2HUt) + H KK (t }] , (11) 

which is obtained [lij by matching the prediction of chiral perturbation theory at 0(p i ) with the contribution of the 
p(770) in the resonance region. The expressions of r^ p (t) (which differs from our decay rate given in section 2) and 
of the loops functions Hpp/(t) can be found in refs. As in our model, the form factor in eq. (11) gives a good 

description of experimental data for the two-pion spectra in the region below Vt = 1.1 GeV. 

The branching ratio for the non-radiative decay that is obtained using the form factor in eq. (11) also underestimates 
the experimental value since: 



B(t -> ttttv) = 21.19% 



(12) 



5 



This low branching ratio reflects again the fact that the form factor in eq. (11) underestimates experimental data of 
the pion form factor for large values of t. One possibility to account for this discrepancy in the predictions of our 
model is to normalize our results for radiative decays in terms of the non-radiative rate. However, for the purposes 
of comparing our results with those of ref. |jj we keep the one-resonance model with the p~ contribution in the 
evaluation of the model-independent radiative amplitudes. 



4. RADIATIVE DECAY MODE 

The Feynman diagrams that contribute to the radiative T~(p) — > it~ (p-)ir (po)v T (q)j(k, e) decay in our meson 
dominance model are shown in Figure J2J). The particles four-momenta are indicated within parenthesis, with k(e) 
denoting the momentum and polarization four-vectors of the photon. 

The decay amplitude has the following generic expansion in powers of the photon energy E 1 [lfij | : 

M = 4- + BE& + CE 1 + • • • , (13) 

where the ellipsis denotes the terms of higher order in E 1 . As we will see below, the terms of order up to E® (Low's 
amplitude) contains only model-independent contributions, while the terms starting at order E 7 arise from model- 
dependent contributions and from the magnetic dipole ((3(0)) and electric quadrupole moments (Q p ) of the p~ meson 
(in this paper we do not consider the possible effects of Q p ). In the following we consider the different contributions 
in more detail. 



A. Model-independent contributions 

The model-independent contributions (Fig. 2.(a-d)) are obtained by attaching the photon to all the charged 
lines and vertices with derivate couplings in the non-radiative Feynman diagram. This set of diagrams leads to a 
gauge-invariant amplitude. In our model, gauge invariance is guaranteed owing to the Ward identity satisfied by the 
electromagnetic coupling and propagator of the p~ introduced in section 2. In other words, we do not need to impose 
gauge invariance to the model-independent amplitude due to the finite width effects of the p~ vector meson. Just for 
a later comparison, let us mention that the effects of the p~-meson magnetic moment (3(0), a gauge-invariant term 
by itself, can not be obtained by imposing gauge invariance to the sum of amplitudes obtained from diagrams (a, c, d) 
in Figure 2. 

Using the Feynman rules corresponding to the vertices and propagators in diagrams (a-d) from Figure 2, we obtain 
the following amplitudes: 



M a = eG F V ud 1 -^^- '"" 2p y —u(qh (l - To )(t>- A + ™ T ) fu(p) (14) 

M b = eG F V ud 9 ^^<j ) --p Q ) u D^(p_+po)T Km D 6 p r L( P -q)e* K l s (15) 

M c = -eG F V ud i9p9 %* P - ' C * (fc + p- - Po ) v D S p l (p-q)k (16) 
V2 p- ■ k p 

M d = eG F V ud ^^e*D 6 o l(p-q)l 5 . (17) 



(P- - Po)kD o -{P- +Pq) 



G F V ud ^^e;D°;_(p- q )l s 



Owing to the Ward Identity given in section 2, it is easy to verify that the model- independent amplitude Mmi = 
M a + Mb + M c + Md is gauge-invariant, namely Mmi{(* — * k) = 0. The amplitude Mmi differs from the 
corresponding model-independent amplitude of ref. |3] in terms of order k and due to the effects of the magnetic 
dipole moment of the p~ meson. As we will see later, the effects of (3(0) are negligible in the integrated observables of 
this radiative decay. However, as it was discussed elsewhere 0:^3' ^s effects can be enhanced with an special choice 
of the kinematics (see section 5.D). 

Just to end this section, we provide the Low's amplitude obtained from Eqs. (14-17) after expanding the amplitude 
Mmi around the soft-photon limit (the form factor f+(t) used here corresponds to the expression in eq. (7) when 



6 




FIG. 2: Feynman diagrams of the model-independent (a-d) and model-dependent (e-k) contributions to r — > -k tx°v t ') decays. 
a = 0): 

M Low = ^^U(t)(^^-^^(p^-po)^ v + ^u(<l)^-- fa) /c f(l- l5 )u(p) 

-/+(*) U u - j^ kV ) l » + {i~ir k ■ po ~ e * ■ po ) {p - p ° riv } ' (18) 

As it can be easily checked, this amplitude coincides with the one obtained in Ref. Q. As is dictated by Low's soft- 
photon theorem [16j, the amplitude depends only on the non-radiative amplitude and on the static electromagnetic 
properties of the external particles. 



B. Model-dependent contributions 

The model-dependent contributions that appear within our meson dominance model are shown in Figures 2. (e-k). 
The diagrams (e-g) contribute to an effective vector hadronic current, while the diagrams (h-k) give rise to an effective 
axial current. We can write these vector and axial model-dependent contributions to the amplitude as follows: 



M A = eG F V ud e* 



A 



fj,6 



4s 



4s 



(19) 
(20) 
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The explicit expressions for the vector and axial terms of the hadronic vertex are the following: 

ye g = gggg^gyo^ [k . pog ^ _ D £ {k + po) [{k + Pq) . {p _ ^ _ {k+ Mp _ ^ 

X D p _*(p-q) (21) 

V f s = g P g pa ^h ai . [k pg ^ kxp _^ D ^ (k +p _ )[{k+p _ ) . {p _ ^ _ {k + p _ Up _ q)k] 

xD^ip-q) (22) 
Vfr = 9p9p ^^ e AWo A V 'D:\k+po)e n , VK , KP y(k+por'D p _ s <\p-q) (23) 

Ks = fjMw e v , m , fl p ^'^'D;'l(k+ Po ){p- q + p^ K (j>- q )s (24) 
V2e[(p — q) A — ' 

A U = f«9 P ™9p*i e n , m ,»pyk»'D^{k+p-){p- q + p Q Up- q ) 5 (25) 
V2e[(p- q ) 2 - m£j p 

^ = '^^ evWMW ,y^(fc +ff ,) [(fc +po ) . (p _ g)g;?B _ (fc +Po ),(p - q) K ] D a y(p - q) (26) 
A% = lfai9p ^ 9p ^ e x , x ,,,p^'k»'D^(k+p-) [(k + p.) ■ (p q)g VK (fc + p-)„(p - q) n ] D ai />(p - q) . (27) 

All the couplings constants appearing in the above expressions can be easily identified from the corresponding 
diagrams in Figure 2. Their values can be fixed from measured decays of the ai, 7r, p and w mesons and will be 
provided in the next section. Note that when the vector and axial vector mesons become heavy degrees of freedom, 
these model-dependent contributions vanish as required by chiral symmetry 

As already anticipated, the above amplitudes are of order one in the photon four-momentum fc. Moreover, they are 
individually gauge-invariant since the conditions k^VH = k^A^g — are satisfied. Of course, the vector and axial 
amplitudes _given above can be decomposed in terms of a basis of four- independent vector and axial tensors as pointed 
out in ref. 



5. DECAY OBSERVABLES 



As it was discussed in the previous section, the decay amplitude of the radiative r lepton decay depends on a 
large set of parameters (coupling constants and masses of mesons) . T he p arameters m p - , g^n entering the model- 
independent contributions where fixed from a fit to experimental data j 1 31 ] in the di-pion mass spectrum of the decay 
r — ► mris in the region below 1.1 GeV. Their values were given in eq. (9) of section 3. The other free parameter 
entering the model-independent amplitude, namely the magnetic dipole moment /3(0) of the p~ meson, is left as a 
free parameter in order to study their effects on the different observables of radiative r decays. 

The model-dependent contributions depend mainly on the values of the coupling constants and masses of the vector 
(p, w), axial a\ and pseudoscalar 7r mesons (we assume isospin symmetry for the masses of neutral and charged 
states) . The values of the coupling constants can be obtained from the measured decay rates of these mesons (except 
for the weak coupling of the a-^ meson whose value is fixed from the Weinberg sum rule f ai = g p ^3)- As ^ was 
the case for the non-radiative decay in section 3, we expect that these effective couplings will give a good estimate 
of the correct size for model-dependent effects. Based on the experimental data compiled in ref. 01 > we wm use the 
following central values: 



g p = 167765.48 MeV (28) 

U = 130.7 MeV (29) 

9pQiT = 4.843 x 1(T 3 MeV -1 (30) 

? 70l7r = 2.9265 x 10~ 4 Me^ 1 (31) 

g pu7T = 0.012 MeV" 1 (32) 

g 7UJ7r = 7.1126 x 10~ 4 Me^ 1 (33) 

g p7ry = 2.2092 x 10~ 4 MeV" 1 . (34) 
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Once we have fixed the values of these parameters, we proceed to compute the different observables associated to 
the radiative r lepton decay. Using the choice of kinematical variables described in the Appendix, we can write the 
differential decay rate in terms of the five independent kinematical variables as follows: 

dr = o ( a\1 2 \ I M. 1 2 dxdtdEy d cos 9 v -d^- , (35) 



where /3_ = ^/T- 4m \ It is the magnitude of the pion velocity in the di-pion rest frame. In order to integrate over 
the relevant kinematical variables we have used the VEGAS [19| integration routine. Next we focus on the results 
obtained for each one of the computed observables. 

A. Branching ratios 

In this subsection we compute the predictions of our model for the branching ratios. As it was done in ref. [jj, 
we distinguish between the bremsstrahlung (model-independent) and the full (that includes also model-dependent 
terms) contributions to the decay rate. Since the unpolarizcd probability is divergent for soft photons, we introduce a 
cutoff energy E™ m to regularize the integral. In Figure 3 we show the branching ratio as a function of E™ m , for the 
normal value of the magnetic dipole moment (3(0) = 2. We compare our results with the predictions based on chiral 
perturbation theory [3j (the three squares in Figure 3). 

We observe that if we exclude the uj meson contribution, diagram in Figure 2.g, we find a very good agreement with 
the calculation of reference Q. However, according to the VMD model, the contribution of the w(782) vector meson 
can not be excluded. This particular model-dependent contribution becomes large due to a particular kinematic 
accident, namely the almost degeneracy of the p~ — uj system, and due to the small decay width of the uj meson 
(T^ = 8.44 MeV). This double resonance effect produces an enhancement of the decay amplitude in approximately 
the same kinematical region. In order to verify this explanation we have increased the mass difference of the p~ uj 
mesons and/or the width of the uj meson and have found that the large effect of the uj meson is decreased in an 
important way. We find that for photon cutoff energies of order E™™ = 200 MeV, the contribution of the uj meson 
becomes already twice the value of all other contributions. Therefore, a measurement of the radiative decay branching 
ratio can help to discriminate among the two models. 

The branching ratio is almost insensitive to reasonable variations in the value of the p~ -meson magnetic dipole 
moment (3(0). In Figure 4 we show the full branching ratio as a function of £™ n for (3(0) = 1,2 and 3. Thus, it is 
clear that this observable cannot help to discriminate values of the magnetic dipole moment. 

B. Photon spectrum 

The photon spectrum can be obtained after integrating over all the kinematical variables in eq. (35) except E 7 . 
This spectrum dT/dE 7 is plotted in Figure 5 for (3(0) = 2. The effect of the uj meson contribution is particularly 
important for E~ > 180 MeV. As in the case of the branching ratio, the photon spectrum is not sensitive to the value 
of the p~ magnetic dipole moment. 

C. Di-pion invariant mass distribution 

Another important observable associated to t~ — > tt^tt^vj is the invariant mass distribution of the di-pion system. 
In the case of the corresponding non-radiative decay, this spectrum shows explicitly the peaks associated to the 
production of vector resonances. It is interesting to study how they are modified by the radiation of photons. On 
another hand, a detailed study of this spectrum in radiative decays is very important in order to remove the hard 
bremsstrahlung from photon inclusive measurements of t~ — > 7r~7r°z/(7) decays |l|. 

In Figure 6 we plot the combined photon and di-pion invariant mass spectra dT /T nr dE 7 dt (T nr is the non-radiative 
decay rate) by choosing (3(0) = 2. In order to avoid the infrared divergences due to the emission of soft-photons, 
we plot our results for a finite value of the photon energy E 7 . Once again, we observe that the presence of the 
w-meson contribution changes the spectrum in a sizable way in all the region of t. However, the position of the peaks 
associated to the photon emission off the ir~ and the t~ external lines are not affected. We also plot in Figure 7 the 
di-pion invariant mass distribution after integrating the previous result for photons of energy larger than 300 MeV 
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FIG. 3: Branching ratio of r — > ir~TT a wy as a function of the soft-photon cutoff E™ n for ,0(0) = 2. The dashed-line denotes 
the model-independent contributions and the solid-line the full contributions. The dotted-line is obtained by excluding the 
contribution of the w meson, diagram 2.g. The points at E-, — 100, 200 and 300 MeV correspond to the full contributions of 
ref. 0). 
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FIG. 4: Branching ratio of r — *• tt ir°wy as a function of the soft-photon cutoff E" 
solid and dotted lines). Only the full contributions are plotted in this case. 



for /3(0) = 1,2,3 (respectively, dashed, 



10 



> 

CD 



1.5 



T3 



0.5 



100 



200 



300 



400 

E y (MeV) 



500 



600 



700 



FIG. 5: Photon spectrum in the r — > -K~-K a v) decay. The solid-line denotes the model-independent contributions, while the 
dotted-line is used for the full contributions. The dashed-line (almost superposed over the solid-line) corresponds to the full 
contribution obtained by excluding fig. 2.g. The observable is normalized to the non-radiative decay rate. 



(unfortunately, a quantitative comparison with results of reference |2( is not possible since they give their results in 
arbitrary units). As in the case of the branching ratios and of the photon spectrum studied in previous subsections, 
the di-pion invariant mass distribution can also help to distinguish between the present model and the one based in 
chiral perturbation theory 



D. Angular and energy photon spectra 

Previous studies of radiative decays involving the production and decay of and on-shell charged vector meson [l7| 
have shown that the angular and energy photon spectra are sensitive to the effects of the vector-meson magnetic 
dipole moment when photons are emitted at small angles. Therefore, we also compute this observable for the case of 
r — > rnrvj decays. 

In Figure 8 we plot the angular and energy photon spectra dT/T nr d cos 8 dEj as a function of the photon energy for 
0=10° and 20° (9 is the angle between the photon and the ir~ in the t~ lepton rest frame) and three different values 
of /3(0). We have subtracted the (well known) contribution arising from the pure bremsstrahlung (terms of order 
k~ 2 in the unpolarized probability) in order to make more visible the effect of /3(0). We observe that there are some 
kinematical regions where the sensitivity to /3(0) is increased and it may eventually help to measure this property. 



6. CONCLUSIONS 



In this paper we have considered the radiative two pion decay of the r lepton. This decay mode was considered 
previously in refs. QUI- The new ingredients of the present study include, (a) an electromagnetic gauge-invariant 
description of the model-independent amplitude including intermediate unstable p~ mesons and, (b) a complete 
calculation of the model-dependent contributions using a meson dominance model. In the framework of the present 
model, all the meson states that are allowed to contribute as intermediate particles were included in the calculation 
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FIG. 6: Distribution in the photon energy and the invariant mass of the di-pion system in r — > ir~ir°i"y decays for = 300 
MeV. The dashed-line (dotted-line) denotes the model-independent (full) contributions. The solid-line is obtained when we 
exclude the meson w, diagram in fig. 2.g. The observable is normalized to the non-radiative rate. 
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FIG. 7: Di-pion invariant mass distribution in r — > ty n°i"y decays for photon energies larger than 300 MeV. Description of 
lines are the same as in Figure 6. The observable is normalized to the non-radiative rate. 



12 




FIG. 8: Reduced angular and energy distribution of photons for three different values of /3(0) = 1,2,3 (dotted, solid and 
dashed lines) and two different angles of the photon emitted with respect to the n~ three-momentum. The model-independent 
contributions (for ,3(0) = 0) have been subtracted. The observable is normalized to the non-radiative rate. 



of the radiative amplitude. In addition, we study the effects of the p~ magnetic dipole moment in the observables of 
the radiative r lepton decay. 

We have found that the branching ratio, the photon spectrum and the di-pion invariant mass spectrum of radiative 
r lepton decays are sensitive to the model-dependent contributions that include an ui meson intermediate state (Figure 
2.g). This contribution produces an important enhancement of these observables with respect to all other contributions 
arising in the present model. The origin of this enhancement can be traced back to the almost degeneracy of the 
p — lo masses and due to the small decay width of the to meson. In the absence of this contribution, our model 
reproduces the results obtained in the framework of chiral perturbation theory Q . Since present formulations of the 
chiral lagrangian interactions do not include the presence of vector-vector-pseudoscalar vertices, it is natural that the 
calculation of reference [jj has not included the contribution of diagram 2.g. Thus, experimental measurements of 
the observables studied in this paper can help to assess the approximation involved in different models. On another 
hand, our calculation confirms that the effects of the model-dependent axial contributions are negligible 0. 

The quantitative difference in model-dependent terms may modify the size of the corrections applied to extract the 
pion form factor from photon inclusive measurements of r — > irirvj decays. As is well known this pion form factor 
seems to be a bit larger than the one obtained from e + e~ — > 7r + 7r _ annihilations for squared momentum transfers 
larger than m 2 pl even after known sources of isospin breaking corrections are taken into account. Since this study 
must consider the effects of virtual radiative corrections, we will consider it elsewhere |2pj . 

The different observables studied in this paper are not sensitive to the effects of the p~ magnetic dipole moment. 
As it was pointed out in ref. 0] the photon emission off the internal charged p meson line (Figure 2.b) is expected 
to carry information about this important property which has not been measured yet. The sensitivity on different 
values of (3(0) is slightly increased when we consider the angular and energy photon spectra for almost collincar 
7r~ — 7 particles 2 . Thus, only processes where the charged p vector- mesons are on their mass-shell can offer a better 



2 As it was concluded in refs. |l7l . choosing such small angles help to suppress radiation off electric charges and make more prominent 
the radiation off the magnetic dipole moment. 
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sensitivity to the magnetic dipole moment [17]. since the radiation off this electromagnetic moment enters at a lower 
order in the photon momentum. 
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APPENDIX: KINEMATICS 

We discuss here the kinematics of the decay T~{p) — > ir~ (p-)ir° (po)v(q)j(k)] for simplicity we choose the isospin 
limit Pq = p_ = m\. The unpolarized squared amplitude for a four-body decay depends upon five independent 
kinematical variables. We choose this set of independent variables to be (we closely follow ref. (21|): 

(x, t, ELy, cos 6^- , fa- ) ■ (36) 

The quantity t = (p_ + po) 2 (x = (q + k) 2 ) denotes the squared invariant mass of the two-pion (1/7) system, _E 7 is 
the photon-energy in the rest frame of the r, and (9^-,4>~) are the spherical coordinate angles that define the tt~ 
three-momentum in the r lepton rest frame. 

The order of the limits of integration can be conveniently chosen according on the energy or angular distribution 
that we want to obtain for the observables. We consider four possible choices (after integrating upon the angular 
variables) : 

• If we integrate successively on _E 7 , t and x, the limits of integration are given by: 

m 2 < x < (m T — 2m 7T ) 2 (37) 

4m^ < t <{m r - ^fx) 2 (38) 

ml + x - t - 2X m 2 T +x-t + 2X 

< L i < zz ( 39 ) 



4m r 4m T 



where 2X = y X(m 2 . 1 x, t), and m 7 is a cutoff parameter introduced to regularize the infrared divergence. 

• If we exchange x <-> t in the order of integration of the previous case, their corresponding limits are: 

Ami < t <{m T - m 7 ) 2 (40) 
m 2 < x <(m T - Vt) 2 . (41) 

• The successive order of integration over t, x and E~,, requires the integration region to be defined as: 

0< x < ^(m;-4m;-2,„ r E T ) 
m T — 2E T 

4m 2 T < t < (m T -2E T )(2m T E 1 -x) (44) 

• Another useful choice is (the order of integration is easily understood) : 

^<_ El <_<^< ,45, 
Ami < 1 <m T (m T -2Ej) (46) 

0< x < 2^-2^-0 
m T — 2Ej 
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Other choices for the order of integrations are also possible. To verify that the different domains of integrations are 
equivalent, we have performed the numerical integrations using the different domains of integration described above 
and have verified that the same results are obtained. 
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